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FUNDAMENTAL PARAMETERS AND CHEMICAL COMPOSITION OF ARCTURUS 
I. Ramirez 1 and C. Allende Prieto 2 ' 3 

To Appear in The Astrophysical Journal 

ABSTRACT 

We derive a self-consistent set of atmospheric parameters and abundances of 17 elements for the 
red giant star Arcturus: T cS = 4286 ± 30 K, logg = 1.66 ± 0.05, and [Fe/H] = -0.52 ± 0.04. The 
effective temperature was determined using model atmosphere fits to the observed spectral energy 
distribution from the blue to the mid- infrared (0.44 to 10 /zm). The surface gravity was calculated 
using the trigonometric parallax of the star and stellar evolution models. A differential abundance 
analysis relative to the solar spectrum allowed us to derive iron abundances from equivalent width 
measurements of 37 Fei and 9 Fell lines, unblended in the spectra of both Arcturus and the Sun; 
the [Fe/H] value adopted is derived from Fei lines. We also determine the mass, radius, and age of 
Arcturus: M = 1.08 ± 0.06 M , R = 25 .4 ± 0.2 i? , and r = 7.1±i;| Gyr. Finally, abundances of the 
following elements are measured from an equivalent width analysis of atomic features: C, O, Na, Mg, 
Al, Si, K, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, and Zn. We find the chemical composition of Arcturus 
typical of that of a local thick-disk star, consistent with its kinematics. 

Subject headings: stars: abundances — stars: fundamental parameters — stars: individual (Arcturus) 



1. INTRODUCTION 

The nearby K-giant Arcturus (HR5340, HD 124897, 
HIP 69673, a Boo) is an excellent reference for spectro- 
scopic studies of giant stars. It is one of the brightest 
stars in the sky (V — —0.05 mag) and its relatively low 
declination (6 ~ +19°) makes it observable from most 
observatories in both hemispheres. Its moderately low 
metallicity ([Fe/H] ~ —0.5) and Galactic space veloci- 
ties associate the star wit h the Milky Way's thick disk 
(e.g.. IRaimrez et al. 2007). Moreover, the star has been 
identified as a member o f a kinematic group. The so- 
called "Arcturus group" (jEggenl fl971[ ) stars have been 
proposed to be members of a dissolved stellar cluster or 
at leas t remnants of a dispersed short-lived star- forming 
event (|De Silva et al.|[2007Hwilliams et al.ll200ll . They 
have been even spec ulated to be of extragalactic origin 
(|Navarro et al.l 120041 ). 

The atmospheric parameters of Arcturus have been es- 
timated by several investigators using a variety of tech- 
niques. The PASTEL database of stellar parameters by 
ISoubiran et al.l (|2010f ). for example, lists 28 entries for 
Arcturus, with 24 of them including [Fe/H] determina- 
tions from high resolution spectra. The simple mean and 
standard deviation for the stellar parameters compiled in 
PASTEL are: T eS = 4324±90K, log.g = 1.71±0.29, and 
[Fe/H] = —0.56 ±0.10. The published parameters do not 
distribute randomly around the mean values due to the 
impact of systematic errors which vary between different 
studies. Therefore, these average literature values are 
not precise. On the other hand, it is possible to obtain 
very precise atmospheric parameters for Arcturus using 
the best quality data and most reliable models available. 
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Moreover, one can attempt to adopt the least model- 
dependent techniques to obtain not only precise but also 
accurate results. 

Having precise and accurate atmospheric parameters 
for Arcturus is important for studies of giant stars in 
general. Since most systematic errors are dependent 
on the atmospheric parameters, differential analyses of 
giant stars relative to Arcturus can largely minimize 
those errors. A dramatic exam ple of the power of dif - 
fere ntial analysis is provi ded by iMelendez et al.l (2009) 
and IRaimrez et al.l (|2009ft . who performed strictly dif- 
ferential analyses of solar twin stars relative to the Sun 
to determine elemental abundances with unprecedented 
precision. The solar parameters arc too different from 
those of the giant stars to be useful for extremely pre- 
cise differential work. Arcturus, on the other hand, 
represents a much better reference for this purpose, 
and some studies have already taken advantage of this 
fact (e.g.. iMcWilliam fc Richl ll994t iWorlev et aD 120091 : 
lAlves-Brito et alJl2010t )~ 

Large spectroscopic surveys are producing enormous 
and homogeneous data sets. By properly analyzing 
them, our knowledge of galaxy formation and evolution 
will be greatly enhanced. We must realize, however, that 
large number statistics does not help removing system- 
atic errors in the spectroscopic analysis, and it is there- 
fore crucial to have reference stars with extremely well 
determined fundamental parameters and abundances. 
The Apache Point Observatory Galactic Evolution Ex- 
periment (APOGEE) survey, in particular, will observe 
about 100,000 giant stars, mainly in the Galactic bulge 



(lAllende Prieto et all 120081: iSchiavon fc Maiewskill201( 
iMajewski et all 120101 : iShetrone et al.l 120101) . A precise 
and accurate determination of fundamental parameters 
and elemental abundances for the giant star Arcturus will 
therefore be of great importance for the proper handling 
of this large data set. 

Arcturus has bee n an alyzed many times in the past (see 
iMackle et aTlll975l and iPeterson et al.lfT993l for just two 
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examples), but the availability of new data, including 
very high signal-to-noise high-resolution spectra, war- 
rants a re-analysis. We use what we consider the most 
reliable methods and models to estimate Arcturus' atmo- 
spheric parameters and elemental abundances. Our work 
is based on classical static 1D-LTE model atmospheres 
and spectrum synthesis, but it sets a good starting point 
for future investi gations on the impact of 3D and non- 



tig_ 

LTE effects (e.g.,|Asriund|[2005) 



2. FUNDAMENTAL PARAMETERS 
2.1. Binarity 

Arcturus is flagg ed in the Hipparcos ca talog as a two 
component object (iPerrvman et al.lll997f) but the adap- 
tive optics observations bv lTurner et al.l (119991) su g gested 
that the star is single. Later. IVerhoelst etaLl (2005) 
found that a binary model with a G-type subgiant sec- 
ondary matched their near-infrared interferometric data 
better than a single star model, implying that the Hip- 
parcos flag is accurate. We find no detectable signatures 
of binarity in the high quality visible and near-infrared 
(normalized) spectrum of the star. 

Arcturus h as a chromosphere (see, e.g., 

lAyres fc Linskyl [1975) and that enhances visibly 
the continuum flux at wavelengths shorter than about 
2000 A. Compared to the Kurucz model atmosphere 
predictions for Arcturus, the star's absolutely ca librated 
fluxes in the ultraviolet, as given bv lAvresI (|2010f ). reveal 
an excess at wavelengths between 2000 and 3000 A that 
may not be related to the chromospheric activity of the 
star. This flux excess could be explained with a binary 
model in which the secondary is a slightly evolved warm 
star, such as that suggested by IVerhoelst et all (|2005[) . 
In this paper we study Arcturus as a single star because 
the impact of a possible secondary is only important 
in the ultraviolet, a wavelength region that we do not 
use and does not affect our analysis. We also note that 
there have been repeated reports in the literature about 
problems to understand the spectrum formed in the 
outermost layers of this star, i n particular CO and H 2 Q 
trans itions (e.g., lAvresI 119861: iRvde et al.l 120021: iTsuiil 
[2009T) . 

2.2. Angular Diameter 

A number of interferometric measurements of the 
angular diameter of Arcturus have been published. 
They are listed in Table [T] along with their sources. 
All reported measurements are consistent within the 
uncertainties. The angular diameter listed in the 
iRichichi et all (f2009h paper, however, has a very small 
error compared to other sources, but it is likely un- 
derestimated, and the same, to less extent, is true for 
the result by iPerrin et al.1 (|1998| ) . This is most likely 
due to the fact that these reported values did not take 
systematics into account and are based on many obser- 
vations, whic h reduces significan t ly the internal errors. 
According to iMozurkewich et all (|2003l ) , systematic un- 
certainties alone limit the precision of these measure- 
ments to no better than 1 %, which is about 0.2 mas 
for Arcturus. Interestingly, this number is very similar 
to the average error bar published by the other groups 
listed in Table pi Therefore, we assu med that the er- 
ror bars for the TRichichi et ail (|2009l ) and IPerrin et al.l 



TABLE 1 
Angular Diameter of Arcturus 

Reference - Facility 



di Benedet to fc Fov (1986) - I2T 
Perrin et al. (1998) - IOTA 
Quirrenbach et al. (1996) - Mklll 
Mozurkewich et al. 7fuu3) - Mklll 
Verhoelst et al"T2005T ^IOTA 
Lacour et aT"72008T ^ IOTA 
20.924 0.003 Richichi et al. (2009) - VLTI 
21.06 0.17 Weighted mean (value adopted) 1 

1 This value was obtained ad opting a more cons erva- 
tive error of . 20 ma s for the IPerrin et al] Q199SI ) and 
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Richichi ct al. (2009) diameters. 



( 1998) measurements are comparable to that from other 
sources, i.e., 0.2 mas. Hereafter, the weighted mean 
and standard deviation from all the values given in Ta- 
ble [1] is adopted as the angular diameter of Arcturus: 
6» LD = 21.06 ±0.17mas. 

Note that the values given in Table[TJ and therefore the 
one adopted here, have been obtained from uniform disk 
measurements combined with limb-darkening corrections 
based on plane-parallel 1D-LTE model atmosphere pre- 
dictions. Our adopted value for the angular diameter is 
thus consistent with the rest of our work, which is en- 
tirely based on the same type of atmospheric models. 

The limb-darkening corrections adopted by each study 
listed in Table [1] are not all from the same source. To 
investigate whether this inconsistency is adding scatter 
to the average value that we adopt, we applied th e limb - 
darkening corrections computed by iDavis et"al"1 (2000) 
using Kurucz model atmospheres to all of the uniform- 
disk inBasyrrenients found in the literature. We used 
the IDavis et al.l corrections as a function of wavelength 
for T eff = 4250 K, logs = 1.5, and [Fe/H] = -0.5, 
which correspond to the node of their grid that is near- 
est to the stellar parameters that we derive for Arc- 
turus. As before, we adopted a minimum error bar 
of 0.2 mas for the published diameters. We obtained 
#ld = 21.04 ± 0.20 mas (see Fig. [1]), which is consistent 
with the value we adopted but has a larger scatter. Thus, 
averaging true angular diameters computed with differ- 
ent limb-darkening corrections does not introduce scatter 
compared to the case when uniform-disk measurements 
from different sources are corrected for limb-darkening in 
a consistent manner. 

Although, based on what can be found in the liter- 
ature, the angular diameter of Arcturus appears to be 
known with high precision (better than 1 % after av- 
eraging independent measurements), we should remark 
that these results are model-dependent. For example, 
the use of spherically symmetric atmospheres instead of 
plane-parallel models results in an incre ase of ~ 0.5 mas 
for th e angular diameter according to IVerhoelst et al.l 
(2005). Moreover, the atmospheres of evolved stars are 
dynamically unstable and their surfaces are not spher- 
ical. This will have an impact on the measurement of 
uniform-disk diameters, which assume symmetric shapes, 
and limb- darkening corrections wh ich are based on model 
atmosphere calculatio ns (see, e.g.. iKoesterke et "al"1l2008l : 
iChiavassa et al.l 120101 ) . In fact, 3D corrections to the 
limb-darkening may compensate the increase of angu- 
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Fig. 1. — Angular diameter measurements as a function of wave- 
length. Uniform-disk measurements are shown with open circles 
(sources are listed in Ta blefTjl. Limb-darken ed diameters, obtained 
using the corrections by Davis ct al. (2000), are shown with filled 
circles. The dashed line shows the weighted average and the dotted 
lines are the ±1 a limits 

lar diameter suggested by static spherical models (e.g., 
lAllende Prieto et al.l 12003 ). This possible systematic er- 
ror will propagate in our analysis but not change our 
results significantly, as shown later. Unless otherwise 
noted, the angular diameter adopted in this work is that 
given in Table [TJ 

2.3. Effective Temperature 

The shape of the spectral energy distribution (SED) 
of a star is determined primarily by its effective temper- 
ature. Certain features are also sensitive to other stel- 
lar parameters such as logy (e.g., the Balmer jump) or 
[Fe/H] (e.g., the G-band and Ca II H&K regions). If 
the latter can be determined accurately using indepen- 
dent methods, however, model fits to the observed SED 
can be used to esti mate T c s with high precision (e.g., 
iRarmrez et~aT1l2006D . 

We use least-squares minimization to estimate the 
effective temperature of Arcturus from model fits to 
spectrophotometric data in the visible, near-infrared, 
and near- to mid-infrared, as described in detail be- 
low. The theoretical fluxes employed are from the Ku- 
rucz grid of no-overshoot model atmospheres with a- 
element enhanced composition ([a/Fe] = +0.4; e.g., 
iCastelli fc Kuruczl [20031) PI The use of plane-parallel 
model atmospheres fo r Arcturus is well justified here 
(jVerhoelst et al.1 120051 ) , although sphericity could have 
an impact on the derivation of the angular diameter from 
interferometric measurements. The observed SEDs are 
divided by the scale factor s = # 2 D /4 so that they rep- 
resent the flux that emerges from the surface of the star. 
For each data set, the spectra (observed or theoretical) 
were smoothed to a common spectral resolution (the one 
that corresponds to the lower resolution spectrum). A 
surface gravity logg = 1.66 (Sect. 12. 5p and iron abun- 
dance [Fe/H] = —0.52 (Sect. were adopted for these 



4 Available online at Ihttp : //kurucz . harvard . edu/gr ids . html I 
These model fluxes have a spectral resolution 01 H = A/Z1A = 
100 — 500 and wavelength sampling with steps of 20 A in the visi- 
ble and 50-200 A at A > 1 fim. 



fits. 

As usual, if the grid of model atmosphere fluxes is given 
by ipx(T c g), the reduced \ 2 value is: 



xl = 



1 



(fx ~ M 2 



n - 2 ^ (A/ A 



(1) 



where fx are the observed fluxes, A /a their errors, and n 
the number of observational data points. If available, we 
used the published A /a values, otherwise we adopted a 
constant 2 % error. The minimization of a set of xt va l~ 
ues allowed us to obtain the best-fit T e ff. We computed 
X 2 U values in steps of 1 K fro m T off = 4220 to 4350 K 
for the visible fits (Sect. I2.3.1[) and in steps of 50 K from 
T cff = 390 to 4600 K for the infrared fits (Sects. 12X21 
and I2.3.3P . To obtain a more precise effective temper- 
ature, we fitted a parabola to the 7 points nearest the 
minimum \t- 

Our least-squares minimization scheme uses only T e ff 
as a free parameter. Surface gravity and iron abundance 
are kept constant. The error in our T e ff estimate is thus 
given by: 



AT eff = 



cir 



\ 1/2 


dT cS 


) + 


ds 



As 



(2) 



where the second term on the right-hand side of this 
equation corresponds to the error introduced by the un- 
certainty in the scale factor s = 6* 2 d/4 and therefore de- 
pends on the precision of the angular diameter measure- 
ment. The errors introduced by the uncertainties in logg 
and [Fe/H] are very small compared to the other sources 
of error. Independently of the input data set, we find 
AT e ff ^ if K using our derived values of A log g = ±0.06 
and A[Fe/H] = ±0.04. Nevertheless, this small error was 
added linearly to each AT e ff estimate. Note that because 
we consider s, logg, and [Fe/H] fixed in the calculation of 
X 2 ,, their contributions to AT c g are added linearly; the 
X 2 , values depend on the adopted values of these fixed 
parameters. 

We note that the Kurucz theoretical fluxes used here 
are not the result of radiative transfer calculations with a 
frequency step fine enough to guarantee that all spectral 
lines are well sampled. Nevertheless, we performed tests 
that show that this will not lead to a significant offset in 
the derived effective temperature. 

2.3.1. Visible 

We us e three publishe d obser ved visible SEDs : iBreeerl 
((1976T) . IKiehlinel (fl987l) . and lAlekseeva etUI (|1996h . 
The Kurucz m odel fluxes have a fine r sam pling than the 
IBreeerl (|1976ft and lAlekseeva et alj (|1996l) data so the 
resolution of the theoretical spectra was degraded when 
obtaining T e a from those d ata sets. The opposite was 
the case for IKiehlingl ( 19871 ). The data and b est model 
fits are shown in Fig. O The IKiehlingl (1987) data are 
severely affected by strong atmospheric absorption fea- 
tures such as O2 and H2O, which are very important 
at wavelengths longer than about 6800 A. We excluded 
the se regions from the mod el fits. At shorter wavelengths 
the lAlekseeva et a"L1 (|1996l) data appear to underestimate 
the UV-blue fluxes. Note, however, that in several other 
regions the fluxes from this source also show important 
discrepancies with the best model predictions. Both the 
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Fig. 2. — Spectral energy distribution of Arcturus, scaled to rep- 
resent the monochromatic flux that emerges from its surface. Open 
circles and solid lines represent observed data (references are given 
in the lower right of each panel). Dotted lines correspond to the 
best fit model atmosphere fluxes. 

iBregerl (| 19761) and lAlekseeva et~aH (|1996f ) fluxes appear 
too high relative to the best fit models near A = 1 fim. To 
ensure consistency in the model fits to these three visible 
data sets, we used a single sample region from 4400 to 
6800 A. The outcome of our calculations does not change 
significantly if wider regions are used, but the model fits 
do not perform as well as in this wavelength window. 

The run of \t values as a function of T g for the three 
observed fluxes is shown in Fig. [21 The best fit T e g values 
for the three observed data sets are consistent within a 
few degrees. Thus, we adopt a weighted average value as 
representative of the T e g that corresponds to the SED in 
the visible: T cB (visible) = 4288 ± 17K. 

2.3.2. Near-IR (1~4 (tm) 

The spectral energy distribution of Arcturus in the 
near infrared, as measured b y NASA's Infrared T elescope 
Facility (IRTF), is given bv lRavner et al.l (|2009D . These 
data cover the region from about 1 to 4 /um, exclud- 
ing two relatively small windows around 1.9 and 2.7 /xm, 
which were discarded due to st rong telluric absorption 
(cf. Fig. 3 in lRavner et a l. 2009). We degraded the spec- 
tral resolution of these data (2000 < R < 2 500) to match 
that o f the Kurucz model fluxes. The iRavner et al.l 
(2009) spectrum has bee n absolutely flux calibr ated us- 
ing 2MASS photometry (jSkrutskie et all 120061) . Note, 
however, that 2MASS photometry is heavily saturated 
in the case of Arcturus and the value listed in the cat- 
alog was obtained by fitting the wings of the point- 
spread- function instead of the full PSF, as was the case 
for most other fainter stars. The uncertainties in the 









o a 




a o & ; 




A r* 


r 






\§> A96 //* 




is \A M 








B76 - 



4250 4300 4350 

Ten (K) 

Fig. 3. — Reduced x 2 value as a function of model T c g for the 
three visible data sets. The exact wavelength range used in these 
computations is 4400-6800 A. Solid, dotted, and dashed lines cor- 
respond to parabolic fits of the 7 points closest to the minima. 
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Fig. 4. — Reduced \ 2 value as a function of model T c g for the 
IRTF data set multiplied by different flux scale factors from 1.00 
to 1.18. The minimum in each case is shown with an open circle. 

2MASS photometry could be underestimated for this ob- 
ject. Thus, we introduced a second free parameter in 
our least-squares fits, namely a constant factor to mul- 
tiply the fluxes, in order to take into account a possi- 
ble large systematic error introduced by the uncertain 
2MASS photometry. 

In Fig.2]we show the xt values as a function of model 
T e ff for flux scale factors from 1.00 to 1.18. Clearly, 
a better fit is obtained after the published fluxes are 
multiplied by 1.12. The best fit model in this case has 
T eff (IRTF) = 4347 ± 69 K. The fact that the published 
fluxes need to be scaled up by 12 % in order to be more 
consistent with the models and the observed visible SED 
suggests that the published 2MASS photometry for Arc- 
turus is indeed uncertain and should be avoided. 

The IRTF fluxes are shown in Fig. [5] for various values 
of the empirical flux scale factor. It is clear that the 
published fluxes (un-scaled) are too low compared to a 
model with the T e g that corresponds to the visible fit 
(Fig. Ek) • The best fit T eS for this case (= 4128 K) is low 
compared to that obtained with the visible SEDs and 
severely underestimates the fluxes there (Fig. [30 . Our 
best fit model (factor=1.12 and T e g — 4347K) does a 
much better job in reproducing the fluxes from about 1.5 
to 4 /xm but overestimates the visible fluxes (Fig. [3)- We 
also note that a factor of 1.08 corresponds to a T e g very 
close to that derived from the visible SEDs. However, the 
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Fig. 5. — Spectral energy distribution of Arcturus from the blue 
to the near infrared. Solid lines c orrespond to t he IRTF data. 
Dotted lines are the visible SED by Brcgcr (1976). Dashed lines 
correspond to the model atmosphere predictions (T e g is given in 
the upper right corner of each panel) . In each panel the IRTF fluxes 
have been scaled by the factor shown in the upper right corner. 

visible and near infrared SEDs do not connect smoothly 
around 1 /im (Fig. [5ji). Excluding the region from 0.8 
to 1.5 /im, a T c ff = 4286 K fits reasonably well both the 
visible and the near infrared data. Hereafter, whenever 
they are used, the IRTF fluxes are multiplied by 1.08. 
However, we stress the fact that the best model fit to 
the IRTF data alone is obtained with a scale factor of 
1.12 and T cff (IRTF) = 4347 ± 69 K, and this is the value 
adopted when averaging out the temperatures derived 
from different spectral windows. 

2.3.3. Near- to Mtd-IR (2-10 fim) 

In this spect ral region, abso l utely calibrated fluxes are 
available from lEngelke et al . (2006). These data were 
smoothed to match the spectral r e soluti on of the theo- 
retical fluxes. The lEngelke et all (|2006l ) SED is based 
on data taken with the Short Wavelength Spectrometer 
(SWS) on the Infrared Space Observatory. Fig. [6] shows 
the SWS data along with visible and near-IR fluxes. The 
theoretical fluxes have a good wavelength coverage up to 
10 /im but it becomes very sparse for longer wavelengths. 
In fact, between 10 and 40 /im, only two theoretical data 
points are available for comparison with the observed 
data. Our least-squares minimization was restricted to 
the range from 2 to 10 /im. Visual inspection reveals that 




1 10 
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Fig. 6. — Spectral energy distribution of Arcturus from t he blue 
to the mid-infrared. Visible data (open circles) are from IBreeerl 
(1976), near infrared data (dashed line) are from IRTF, and the 
near- to mid-infrared data (thick solid line) are from SWS. The 
dotted line corresponds to the best fit model. The bottom left 
corner inset shows an expanded view of the 0.6 to 2.0 (im region. 

the two model flux points beyond 10 /im are fully con- 
sistent with the best model obtained from the shorter 
wavelength data. The best fit model was obtained, as 
usual, with the least-squares technique, and found to be 
T eff (SWS) = 4152 ± 84 K. This effective temperature is 
significantly cooler than that obtained with the shorter 
wavelength data. Even though the error of this value 
is the largest of the three T e g estimates, it is clear that 
the best fit model to the SWS data underestimates the 
visible and near-IR data, i.e., the fluxes at wavelengths 
shorter than about 2 /im, as shown in the bottom left 
corner inset of Fig. |6] 

2.3.4. Adopted T cB 

The effective temperature of Arcturus, as suggested 
by model fits to observed SEDs in three different spec- 
tral regions, is given in Table [21 The weighted mean and 
standard deviation from these three results is: T c g = 
4286 ± 30 K. In Fig. [7] we show the spectral energy distri- 
bution of Arcturus from near-UV to mid-infrared wave- 
lengths. The observed data are from different sources. In 
the visible, we computed an average of all observations 
availa ble, interpolat ing the data sets to the wavelengths 
of the lBregerl (| 19761 ) data, excluding the iKiehlingl (|1987f ) 
fluxes for wavelength regions affected by telluric absorp- 
tion. A Kurucz' model flux distribution of T e g = 4286 K 
is also shown in Fig. [71 along with similar models for 
r Gff ± 100 K and T cS ± 200 K. Clearly, the sensitivity to 
T e ff is more important at shorter wavelengths and this 
is why the model fits in the near- to mid-infrared re- 
gions have a larger uncertainty and therefore much lower 
weight in the calculation of the final T c ff ■ This figure also 
suggests that our error estimate is realistic. The observed 
fluxes are fully contained within the ±100 K model flux 
limits while inspection of the ±50 K limits shows that 
most of the observed visible data are contained within 
those limits, as shown in the bottom left corner inset of 
Fig. [7] Direct integration of the composite SED shown in 
Fig. [7] results in a bolometric flux (on Arcturus' surface) 
/boi = / f\d\ — 1.923 x 10 10 erg cm -2 s" 1 , which corre- 
sponds to an effective temperature T G ff = (/boi/c) 1 ^ 4 = 
4291 K (er is the Stefan-Boltzmann constant), which is 
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TABLE 2 

Effective Temperature of Arcturus from SED Fits 



Spectral 


T c ft 


error 


region 


(K) 


(K) 


Visible 


4288 


17 


IRTF 


4347 


69 


sws 


4152 


84 


Adopted 


4286 


30 



TABLE 3 

Trigonometric Parallax of Arcturus 



7T 


error 


Reference 


(mas) 


(mas) 




88.5 


2.1 


Harrington ct al. (1993) 


88.4 


1.8 


van Altcna ct al. (1995) 


87.53 


1.48 


Gatewood (2008) 


88.83 


0.54 


van Leeuwen (2007) 


88.65 


0.40 


Weighted mean (value adopted) 



in excellent agreement with our adopted T c g from the 
model fits. 

As explained in Sect. 12.21 our adopted angular diam- 
eter of Arcturus could be affected by a systematic er- 
ror of up to 0.5 mas, which will have an impact on the 
calculations presented in this Section. We repeated the 
SED fits using a larger angular diameter (21.56 mas) and 
found T ff ~ 4260 K, i.e., an effective temperature cooler 
by about 30 K0 We also found slightly larger minimum 
reduced x 2 values in all fits. Considering that a 0.5 mas 
correction is probably an extreme case, we conclude that 
the impact of systematic errors on the angular diameter 
is not crucial for the rest of our analysis. 

2.4. Parallax 

A measurement of the trigonometric parallax of Arc- 
turus is available in the Hi pparcos catalog. T he new 
reduction of Hipparcos data ()van Leeuwenll2007l ) did not 
change significantly the mean value of the parallax (it 
went from 88.85 to 88.83 mas; a distance of 11.3 pc) but 
the error bar decreased by almost 30%. The Hipparcos 
parallaxes are in general much more reliable and precise 
than ground-based measurements. However, Arcturus is 
one of the few very bright stars in the catalog and thus 
the results might suffer from systematic errors. To en- 
sure that the Hipparcos parallax of Arcturus is reliable, 
we checked ground-based measurements. 

Table [3] lists independent measurements of the trigono- 
metric parallax of Arcturus. The first three en- 
trie s are ground-ba s ed. Moreover, the value given 
by Ivan Altena et al] (Qj)95) represents the average of 
6 previously publishe d measurements (not including 
iHarrington et al] Il993j ) , which were all in good agree- 
ment. Although the average of ground-based measure- 
ments appears to be slightly lower (~ 88.1 mas) com- 
pared to the Hipparcos value (88.83 mas), all measure- 
ments are consistent with each other within the esti- 
mated 1 ct errors. This suggests that there are no real 

5 Note that if T e ff were derived directly from the bolometric flux 
and angular diameter, the change in T c g would be larger, about 
50 K. The value derived here (~ 30 K) corresponds to the SED fits, 
which depend not only on the scale factor s but also on the shape 
of the flux distributions, particularly in the visible, the spectral 
region that has largest weight when computing the final T fj. 



problems with the Hipparcos parallaxes of very bright 
stars. Thus, we adopted the weighted average of the 
parallaxes listed in Table [3] as the parallax of Arcturus: 
tt = 88.65 ±0.40 mas. 

2.5. Surface Gravity, Age, and Mass 

Given that we know the effective temperature of Arc- 
turus with high precision (0.7%) and that its parallax is 
accurately known (0.5 %), the most reliable way to deter- 
mine the star's surface gravity (log g) consists on placing 
the star on the HR diagram and comparing its location 
with theoretical predictions based on stellar evolution 
calculations. This approach also allows us to estimate 
the star's mass (M) and age (r). 

We use a classical isochrone fitting technique to de- 
termine Arcturus' age, mass, and log g. Details of our 
particular implementation will be given in a forthcom- 
ing p aper (iRamfrez et al.ll201lL but see also lReddv et al] 
120031 and [Allende Prieto et al]|2004l who use essentially 
the same approa ch). Briefly, we use a fine grid of Yonsei- 
Yale isochrones (|Yi et al J 120011 iKim et all 120021 ). which 
consider a-element enhancement at sub-solar metallici- 
ties, increasing linearly from [a/Fe] = 0.0 at [Fe/H] = 0.0 
to [a/Ye] = +0.3 at [Fe/H] = -1.0. We calculate the 
absolute magnitude of Arcturus from its observed ap- 
parent magnitude and trigonometric parallax. The ap- 
parent magnitude adopted (V — —0.051 ± 0.013) is the 
on e listed in the Gen e ral C atalogue of Photometric Data 
bv iMermilliod et all (|1997l GCPD), a value that corre- 
sponds to a weighted average of more than 89 measure- 
ments from 14 independent sources. We checked this 
value by direct integration of the observed spectral en- 
ergy distributions used in Sec t. 12.3.11 conv o lved w ith the 
V filter response function bv lCohen et al.l (|2003| ). which 
ensures a result on the standard Landolt scale. We ob- 
tained V = -0.049 ± 0.011, i.e., fully consistent with 
the GCPD. For the calculation of the absolute magni- 
tude, apparent magnitude and parallax errors were prop- 
agated linearly. We obtain M v = -0.313 ± 0.016. We 
calculate a probability distribution function (PDF) for 
the isochrone points within 3 a from the observables My, 
T c s, and [Fc/H]. PDFs are computed for the age, mass, 
and surface gravity, as shown in Fig. [8j and the values 
at their maximum are adopted as the stellar parameters. 
The PDFs are generally asymmetric. Therefore, lower 
and upper 1 a and 2 a Gaussian-like error bars can be 
determined from the shape of the PDFs. We adopt the 
1 a values as our limits for the error bar estimates. Given 
the high precision of our derived atmospheric parameters, 
these PDFs have well defined peaks. The parameters 
derived from these distributions are: r = 7.1^ 2 Gyr, 

M = 1.08 ± 0.06 M , and logg = 1.66±g;|. Hereafter we 
adopt log g — 1.66±0.05. We note that these values agree 
fairly well with some of the early spectroscop ic analysis 
in the literature (e.g., lAvres fc Johnsonlll977l and refer- 
ences therein). 

A number of recent studies have pointed out that the 
determination of ages and other stellar parameters from 
isochrones ca n be severely affected by a number of sta- 
tistical biases (iNordstrom et al.ll2004l:lPont fc Everll2004t 
I J0rgensen fc Lindegrenll2005l: Ida Silva et alj|200ffl which, 
can be addressed, for example, usi ng Bayesian methods . 
We did not take this approach but Ida Silva et al.l (2006) 
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Fig. 7, — Near-UV t o mid-infrared spectral energy distribution of Arcturus. Observed data are from Kichling ( 1987), Brcgcr ( 1976), and 
Aleksccva ct al. (1996) for the visible (average values shown), IRTF (solid line from 1.1 to 4/xm), and SWS (solid line from 4 /mi). The 
best fit model is shown with the dotted line. Model fluxes corresponding to T c g ± 100 K and T c g ± 200 K are also shown with dotted lines. 
The bottom left corner inset shows an expanded view of the 0.3 to 1.2 (im region. 
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Fig. 8. — Age, mass, and logg probability distribution functions 
from our isochrone fit. The vertical solid line is located at the 
maximum of the PDFs while the location of the Gaussian-like 1 a 
and 2 cr lower and upper limits are shown with dotted and dashed 
lines, respectively. 



provide an online tool to derive the stellar age, mass, 
and surface gravity using a Bayesian approach from any 



input observational data@ We used this tool and found, 
for our derived parameters of Arcturus: r = 8.4 ± 
2.3 Gyr, M = 1.00 ± 0.09 M Q , and log 5 = 1.67 ± 0.06. 
These values are all consistent with our estimates within 
the quoted 1 a errors. Note that in addition to im- 
proving the met hod of stellar parame ter determination 
from isochrones, Ida Silva et al.l (|2006l ) used a different 
set of isochrones (those from the Pad ova group; e.g., 
iBertelli et aT1ll994l : iGhardi et al.ll2000f l. which could be 
another source of systematic uncertainty. The excellent 
agreement between t he logg values that w e derive and 
those obtained using Ida Silva et ail (|2006l ) implementa- 
tion suggests that the impact of the Bayesian approach 
and the choice of isochrones is relatively small regarding 
the determination of surface gravity (and probably mass 
and age) of stars like Arcturus. 

Using our derived surface gravity and mass from 
isochrones we can calculate the radius of Arcturus: R = 
25.2 ± 0.4 Rq. Alternatively, one could use the direct 
and independent measurements of angular diameter and 
parallax to calculate this radius: R — 25.5 ± 0.2 Rq. 
The good agreement between these two estimates of Arc- 
turus' radius gives us confidence that our derived log g 
value is not only precise but also accurate. Averag- 
ing the two estimates of Arcturus' radius we obtain: 
R = 25 .4 ± 0.2 Rq. If we use an angular diameter larger 
by 0.5 mas (Sect. E2]), we obtain R = 25.6 ± 0.4 R Q and 
R = 26.1±0.2 Rq from the two methods described above. 
The mass and log g value were re-computed consistently 
owing to the lower T e g that corresponds to the larger 
angular diameter. Thus, in this case we obtain slightly 
larger radii, which are, however, not as consistent with 
each other as before. 

2.6. Iron Abundance 

The very high signal-to-noise, high reso lution spectrum 
of Arcturus by Hinkle & Wallace (2005) was inspected 
alongside the even higher signal-to-noise, higher resolu- 

6 http:/ /stev. oapd.inaf.it/cgi-bin/param 
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tion solar spectrum by iKurucz et al.l (|1984l ) in order to 
find spectral lines due to iron (both neutral and singly 
ionized) which are unblended in both spectra and have a 
well defined local continuum. Furthermore, we selected 
only lines that have reliable log gf values measured in the 
laborato ry. Our adopted log o / values have been com- 
piled by iRarmrez et al.l (|201lh . We refer the reader to 
this paper for details on the sources of laboratory atomic 
data. The resulting linelist consists of 37 Fei lines and 
9 Fe II lines, which are given in Table |4j along with the 
relevant atomic data and equivalent widths (EWs) mea- 
sured in the solar and Arcturus' spectra. The latter 
were measured using IRAF's task splotQ fitting Gaus- 
sian profiles to most lines but Voigt profiles to features 
with extended wings. Although both the spectra of the 
Sun and Arcturus are already continuum normalized, in 
each of our EW measurements we used a local pseudo- 
continuum determined by visual inspection of a relatively 
wide spectral window (typically ±5 A around each fea- 
ture). While for the Sun this makes little difference, it 
can be important for certain spectral regions in the spec- 
trum of Arcturus due to the presence of a large number 
of weak molecular features. Using the pseudo-continuum 
minimizes the impact of these small blending features. 

The latest versions of Kurucz' no-overshoot model 
atmos pheres a nd the spectrum synthesis code MOOG 
(e.g., ISnedenl H97l were used to compute elemen- 
tal abundances. For Arcturus, we adopted a model 
atmosphere with a-element enhancement ([a/Fe] = 
+0.4). van der Waals damping cons t ants were 
adopted from the works by IBarklem et al.l (j2000f ) and 
IBarklem fc Aspelund- Johansson! (|2005[) . 

Abundances were measured differentially, on a line-by- 
line basis, with respect to the Sun. In every abundance 
calculation described below, a microturbulent velocity 
(vt) was derived so that the abundances do not correlate 
with the reduced equivalent widths (REW = log EW/X) 
of the Fei lines. The microturbulent velocity is thus dif- 
ferent for each case tested but this has a minor impact 
on the results and the discussion below. 

Using our derived values of T e g — 4286 K and \ogg = 
1.66, we obtain a mean [Fe/H] = —0.52 ± 0.02 from 
the Fei lines and [Fe/H] = -0.40 ± 0.03 from the Fen 
lines (see Fig. [9j the derived microturbulent velocity is 
v t = 1.74 km s _1 ). The error bars here correspond 
only to the line-to-line scatter and do not include the 
uncertainty introduced by errors in T a g and logg. The 
mean iron abundances that we obtain from Fe I and Fe n 
lines separately are thus inconsistent by 0.12 dex. On 
the other hand, the iron abundances inferred from the 
Fei lines do not correlate significantly with excitation 
potential. 

In Fig. [TU] we show the iron abundance inferred from 
Fe I and Fe II lines as a function of T e g for various values 
of log g. Clearly, the Fe n lines are more sensitive to both 
parameters. It is not possible to reconcile the Fei and 
Fell abundances within our 1D-LTE approach for the ef- 
fective temperature derived in Section 12.31 Further work 
on the effects of non-LTE and possibly also surface inho- 

7 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, under cooperative agreement with the Na- 
tional Science Foundation. 



TABLE 4 
Iron Line List 



Wavelength 


EP 


logs/ 


EW Sun 


EW Arcturus 


(A) 


(eV) 




(mA) 


(mA) 


Fei 


5295.3101 


4.420 


-1.590 


29.0 


47.7 


5379.5698 


3.690 


-1.510 


62.5 


99.0 


5386.3301 


4.150 


-1.670 


32.6 


56.0 


5441.3398 


4.310 


-1.630 


32.5 


55.7 


5638.2598 


4.220 


-0.770 


80.0 


104.9 


5679.0229 


4.652 


-0.750 


59.6 


72.7 


5705.4639 


4.301 


-1.355 


38.0 


62.7 


5731.7598 


4.260 


-1.200 


57.7 


83.0 


5778.4531 


2.588 


-3.440 


22.1 


74.7 


5793.9141 


4.220 


-1.619 


34.2 


56.6 


5855.0762 


4.608 


-1.478 


22.4 


37.1 


5905.6699 


4.650 


-0.690 


58.6 


74.7 


5927.7900 


4.650 


-0.990 


42.9 


56.3 


5929.6802 


4.550 


-1.310 


40.0 


57.1 


6003.0098 


3.880 


-1.060 


84.0 


112.4 


6027.0498 


4.076 


-1.090 


64.2 


93.6 


6056.0000 


4.730 


-0.400 


72.6 


85.2 


6079.0098 


4.650 


-1.020 


45.6 


60.2 


6093.6440 


4.607 


-1.300 


30.9 


45.9 


6096.6650 


3.984 


-1.810 


37.6 


64.7 


6151.6182 


2.176 


-3.282 


49.8 


119.7 


6165.3599 


4.143 


-1.460 


44.8 


70.6 


6187.9902 


3.940 


-1.620 


47.6 


78.4 


6240.6460 


2.223 


-3.287 


48.2 


118.7 


6270.2251 


2.858 


-2.540 


52.4 


107.9 


6703.5669 


2.759 


-3.023 


36.8 


92.7 


6705.1021 


4.607 


-0.980 


46.4 


64.3 


6713.7451 


4.795 


-1.400 


21.2 


30.6 


6726.6670 


4.607 


-1.030 


46.9 


62.4 


6793.2588 


4.076 


-2.326 


12.8 


29.4 


6810.2632 


4.607 


-0.986 


50.0 


66.4 


6828.5898 


4.640 


-0.820 


55.9 


72.2 


6842.6899 


4.640 


-1.220 


39.1 


55.0 


6843.6602 


4.550 


-0.830 


60.9 


80.0 


6999.8799 


4.100 


-1.460 


53.9 


78.8 


7022.9502 


4.190 


-1.150 


64.5 


89.9 


7132.9902 


4.080 


-1.650 


43.1 


68.6 


Fen 


4576.3330 


2.844 


-2.950 


64.6 


76.7 


4620.5132 


2.828 


-3.210 


50.4 


60.1 


5234.6240 


3.221 


-2.180 


82.9 


88.7 


5264.8042 


3.230 


-3.130 


46.1 


47.9 


5414.0718 


3.221 


-3.580 


27.3 


31.0 


5425.2568 


3.200 


-3.220 


41.9 


45.4 


6369.4619 


2.891 


-4.110 


19.2 


23.1 


6432.6758 


2.891 


-3.570 


41.3 


47.0 


6516.0771 


2.891 


-3.310 


54.7 


59.0 



mogeneities is needed to tackle this problem, but this is 
beyond the scope of this paper. Ionization balance would 
be achieved if the effective temperature of Arcturus were 
about 4380 K, which is inconsistent with our preferred 
value for this parameter (our derived T e g would have to 
be off by more than 3 a) . 

The loci of stellar parameters T e g and log g for which 
the conditions of ionization and excitation balance are 
satisfied are shown in Fig. QTJ This figure shows that the 
iron abundance analysis with our preferred T c g and log g 
values satisfies excitation balance but not ionization bal- 
ance. Thus, at least within our 1D-LTE approach, the 
analysis of Fe I lines seems more reliable than that for the 
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Fig. 9. — Iron abundance of Arcturus, differential on a line-by-line 
basis with respect to the Sun, as a function of excitation potential 
and reduced equivalent width of the lines. Neutral iron lines (Fel) 
are represented with open circles while singly ionized iron lines 
(Fe II) are shown with filled squares. 
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Fig. 10. — Iron abundance derived from Fel (open squares con- 
nected by solid lines) and Fell (filled circles connected by dashed 
lines) spectral lines as a function of input T c ff. Three possible 
choices of log g = 1.60, 1.66, 1.72 are tested. Thick lines correspond 
to our preferred value of logg. 

Fell lines. Note that at the temperatures found in the 
line-forming layers of photospheres of stars like Arcturus, 
iron is found mostly in its neutral stage, thus making Fe I 
lines more robust against departures from LTE and errors 
in the atmospheric parameters. On the other hand, Fe n 
dominates in the layers where the continuum is formed. 
We adopt the iron abundance from Fe I lines as the iron 
abundance of Arcturus because it is internally more ro- 
bust and consistent with the rest of our 1D-LTE analysis. 
The discrepancy regarding the Fen line analysis cannot 
be solved within the 1D-LTE approach; it must be ad- 
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Fig. 11. — The location of T c ff,logg pairs for which ionization 
balance is satisfied is shown with the solid line. Similarly, the 
location of T e g , logg pairs for which excitation balance is satisfied 
is shown with the dashed line. The filled circle with error bars 
corresponds to our derived values of T e ff and logg for Arcturus, 
obtained independently from the iron line analysis. 

dressed in future work but for now it should be accepted 
as one of the limitations of the standard spectroscopic 
approach for the determination of iron abundances^ 

Hereafter, the value of iron abundance that we adopt 
is [Fe/H] = —0.52 ± 0.04. In this case we have also in- 
cluded the error introduced by the uncertain T e ff and 
log<? values, in addition to the line-by-line scatter. 

The results described above do not change if we use 
absolute abundances instead of differential ones. How- 
ever, the internal errors are reduced significantly if we 
use differential analysis. This suggests that the only ad- 
vantage of using a solar spectrum as reference for dif- 
ferential analysis of Arcturus (and probably other red 
giant stars) is to minimize the errors due to uncertain 
log g/ values. Systematic errors due to simplifications 
in the model atmosphere computations and/or spectral 
line synthesis are not removed with a solar differential 
analysis. 

2.7. Kurucz vs. MARCS Model Atmospheres 

For consistency, in this paper we use only plane- 
parallel Kurucz model atmospheres. However, we re- 
peated the entire procedure described in this section 
to derive the atmospheric paramete rs using the lat- 
est M ARCS model atmosphere grid (jGustafsso n et al.1 
l2008f ) with standard chemical composition, which means 
that a-element enhancement is considered for stars with 
[Fe/H] < 00 Plane-parallel MARCS models are only 

8 The p ossib le systematic error of 0.5 mas in the angular diam- 
eter (Sect. 12. 2|l would not solve these problems. The decrease of 
about 30 K in T e f>f worsens the ionization imbalance and it is not 
compensated enough by the small lowering of the log g value, which 
is about 0.04 dex. 

9 The MARCS grid of atmospheric models used in this work is 
available online at http://marcs.astro.uu.se 
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available for log g > 3.0 so for this analysis of Arc- 
turus we used spherically symmetric MARCS models. 
We find: T eS = 4282 ± 36 K, logg = 1.66 ± 0.05, and 
[Fe/H] = — 0.54± 0.05, i.e., fully consistent with the re- 
sults obtained with Kurucz model atmospheres. The iron 
abundance analysis was also made differentially on a line- 
by-line basis with respect to the Sun. 

The discussion regarding ionization balance using Ku- 
rucz models, given in the previous section, does not 
change qualitatively for MARCS models. The mean 
Fell abundance in this case is also about 0.1 dex higher 
than that from the Fei lines. The only important dif- 
ference with respect to the Kurucz model atmosphere 
analysis is that the excitation balance is no longer sat- 
isfied (hence the slightly large error for [Fe/H]). For a 
log g = 1.66, the effective temperature has to increase to 
about 4370 K if MARCS models are used. Kurucz mod- 
els, on the other hand, are consistent with excitation 
balance if T c s ~ 4290 K, as suggested by the spectral 
energy distribution. Despite this discrepancy, the at- 
mospheric parameters derived are nearly independent on 
which set of the most commonly used model atmosphere 
grids, namely Kurucz and MARCS, are used. Note that 
the excitation balance is heavily sensitive to the 5 Fei 
lines with EP < 3 eV whereas many more high excita- 
tion potential lines are available. Small errors in the 
measured equivalent widths could affect somewhat these 
results. Thus, we must not imply a superiority of one set 
of models over the other based on our Fe I line analysis. 

3. ELEMENTAL ABUNDANCES 

We employed a curve-of-growth (COG) approach, us- 
ing MOOG, to measure elemental abundances of 16 ele- 
ments other than iron. Equivalent widths were measured 
carefully with IRAF's task splot, de-blending lines when 
necessary. For each element other than K, more than 
one line was available. The weighted average abundance 
was finally adopted, with the weights being determined 
from the actual COGs (i.e., from the slope of the abun- 
dance versus reduced equivalent width relation, assum- 
ing a typical error of 1% for the latteiQ). Thus, strong 
and possibly saturated lines, in particular those in the 
spectrum of Arcturus, are given lower weight (in this 
way we also minimize the impact of uncertain broaden- 
ing parameters; see below). Similar to the case of iron, 
we calculated abundances on a line-by-line basis and dif- 
ferentially with respect to the Sun before averaging, in 
order to reduce the impact of errors in the atomic data 
(with a few exceptions, as explained below). This im- 
plies that not all available lines were used but only those 
for which reliable equivalent widths could be measured in 
the spectra of both Arcturus and the Sun. Our measured 
equivalent widths, atomic line data, and absolute abun- 
dances are given in Table [5j The abundance errors in 
this table correspond only to those from the COG anal- 

10 Formal EW errors c an be compute d from the properties of 
the observed spectra, as injpayrcl (1988). We find errors of order 
0.1 mA, which is about 0.2 % for an EW = 50 mA line. We prefer to 
adopt a larger error because of the uncertainties introduced by the 
continuum placement, which roughly scale with line strength. At 
this point the error bars are used only to weight the abundances 
inferred from different lines of the same element so our results 
will not be significantly affected if these errors are under- or over- 
estimated, as long as their relative values are correct. 



ysis and do not include the uncertainties from the stellar 
parameters. 

Inspection of Table [5] reveals that the spectral lines in 
the spectrum of Arcturus tend to be stronger than those 
in the solar spectrum. This is also the case of the iron 
lines employed. Pressure broadening could affect signif- 
icantly the EWs of strong lines while being less impor- 
tant for weak features. Thus, it is necessary to investi- 
gate possible errors introduced by this ingredient on the 
line formation calculations. We re-computed the abun- 
dances using the classical formula by Unsold to calculate 
the van der Waals damping c onstants instead of using 
those by iBarklem et all (|2000| ). On a line- by- line basis, 
the differential iron abundances increased between 0.00 
and 0.04 dex, depending on the line, while the average 
[Fc/H] increased by 0.02 dex. For the abundance ratios, 
[X/Fe], we obtained shifts of only ±0.02 dex. It has been 
shown that the Unsold approximation severely under- 
est imates the damping constants while the calculations 
bv IBarklem et~aH (2000) improve abundance determina- 
tions. Thus, the differences quoted here, albeit small, 
should still not be considered as potential systematic er- 
rors. They only represent the worst case scenario of er- 
rors introduced by uncertain damping constants. Other 
pressure broadening prescriptions will result in abun- 
dance ratios much m ore similar to those obtained using 
IBarklem et all pOOOh constants. 

Only atomic lines have been used for these calculations 
even though the spectra of cool giant stars like Arcturus 
are rich in molecular features. We prefer to avoid the 
latter in our work because of the potentially severe er- 
rors that could be introduced by surface inhomogeneities 
which are not taken into account in the modeling of 
the star's atmosphere (e.g.. I Asplundl [20051 ) . It has been 
shown that these so-called 3D effects are very important 
for molecular features in cool giants, an e ffect that is 
further enh anced by a low metallicity (e.g., Co llet et al.l 
120071 12009ft . While these studies are mostly theoretical, 
some fundamental predictions of these 3D models have 
been tested against high quality observations, showing in 
general good agreement and therefore providing support 
to the findings of la rge 3D abundance corrections (e.g., 
iRami'rez et aJLl[2010t ). 

Hyperfine splitting was taken into account for the anal- 
ysis of Mn and Co lines. The equivalent widths of these 
features could be measured with high precision, so line 
profile fitting was not necessary. Instead, we used the 
"blends" driver in MOOG to take the effect into ac- 
count. Hyperfine structure constants were adopted from 
Kurucf3 but with the transiti on probabilities scale d to 
the total logg/ values given in lAsplund et a l. (2009). 

Nearly all spectra l lines used in t his w ork come from 
the line selection bv lAsplund et all (|2009l ). therefore en- 
suring that our solar reference analysis is as accurate as 
possible. Indeed, the aver age difference b etween our solar 
abundances and those bv lAsplund et all excluding oxy- 
gen (see below), is A(A X ) = -0.014± 0.082. The largest 
differenc e is seen for K , for which we find = 5.31 
whereas lAsplund et al.l derive Ak — 5.03. This element 
has only one very strong line available for analysis so 
its abundance as reported in t his work should not be 
considered very reliable. In fact, iZhang et al.l (|2006l ) re- 

11 http://kurucz.harvard.edu 
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TABLE 5 

Line List for Elements other than Iron 
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8.37 


0.02 


29.0 


8.32 


0.02 


-0.05 


0.03 


C i 


9078.28 


7.48 


-0.57 


106.3 


8.47 


0.02 


33.1 


8.30 


0.02 


-0.17 


0.03 


Ci 


9111.80 


7.49 


-0.30 


128.1 


8.40 


0.02 


45.0 


8.32 


0.02 


-0.08 


0.03 


Oi 


5577.34 


1.97 


-8.24 




8.69 




10.5 


8.68 


0.01 


-0.01 


0.01 


Oi 


6300.30 


0.00 


-9.72 




8.69 




68.1 


8.66 


0.02 


-0.03 


0.02 


Oi 


6363.78 


0.02 


-10.19 




8.69 




31.1 


8.66 


0.01 


-0.03 


0.01 


Nai 


4751.82 


2.10 


-2.08 


11.5 


6.20 


0.01 


29.0 


5.77 


0.01 


-0.43 


0.01 


Nai 


5148.84 


2.10 


-2.04 


12.9 


6.21 


0.01 


35.5 


5.83 


0.01 


-0.38 


0.02 


Nai 


6154.23 


2.10 


-1.55 


37.0 


6.26 


0.01 


73.3 


5.84 


0.02 
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0.03 


Nai 


6160.75 


2.10 


-1.25 


55.8 


6.25 


0.02 


92.9 


5.82 


0.03 


-0.43 


0.03 



port a non-LTE correction of —0.29 for this line, which 
would bring our K abund ance into excell ent agreement 
with the value derived by lAsplund et al.l E xcluding K, 
the average difference with lAsplund et al.l solar abun- 
dances is A(Ax) = —0.033 ± 0.033; i.e., our solar abun- 
dances are marginally lower. We note that, in general, 
th e differences be tween our solar abundances and those 
by lAsplund et al.l can be reconciled if we take non-LTE 
effects into account, as they did 13 The latter are domi- 
nated by over-ionization if weak and moderately strong 
lines are used, resulting in lower abundances if derived 
from neutral species, as is our case (with few exceptions). 

Our adopted atomic data are also from the compila- 
tion bv lAsplund et al.l (|2009). Their line selection is very 
strict, therefore leaving us with some elements with very 
few or no measurable lines available in both the solar and 
Arcturus' spectra. This was the case of C and Al. For 
these elements we complemented the linelist with atomic 
data from the NISTl3 dat abase (for the 8335, 90 78, and 
9112 A C I lines) and from lMelendez et al.l (|2009( ) for the 
7836 A All line. The oxygen lines used here are not 
easy to measure in the solar spectrum because they are 
very weak and heavily blended; detailed line synthesis 
must be employed for those features. For these lines, we 
adopted the a verage solar oxygen abundance given by 
lAsplund et "all (A = 8.69). 

The weighted average abundances relative to the solar 
abundances we derive for Arcturus are listed in Table [6l 
where we also summarize our error analysis. The con- 
tribution to the error by the line-to-line scatter is given 
by (Tiino, and corresponds to the standard error of the 
line-by-line abundance dispersion for each element. The 
errors introduced by our formal uncertainties in the at- 
mospheric parameters are given by at (corresponding to 
the T e ff uncertainty), a g (for logg), and a m (for [Fe/H]). 
They were computed using the abundances derived from 
model atmospheres with slightly modified stellar param- 
eters compared to those obtained with our preferred val- 
ues. The four error contributions were then added in 
quadrature to obtain the final error estimates, which are 
given in the last column of Table [6] This is equivalent to 
assuming that the covariances between the uncertainties 
in the atmospheric parameters are small, which is ap- 

12 They also consider 3D effects but the 3D-1D differences are 
modest for atomic lines. 

13 http://physics.nist.gov/PhysRefData/ASD/lines_form.html 



TABLE 6 

Error Analysis and Mean Elemental Abundances of 
Arcturus 
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Ci 


-0.09 


0.03 
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0.035 


0.005 


0.43 ± 0.07 


Oi 


-0.02 


0.02 
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0.020 


0.010 


0.50 ±0.03 


Nai 


-0.41 


0.02 


0.025 


0.000 


0.000 


0.11 ± 0.03 


Mgi 


-0.15 


0.03 


0.005 


0.010 


0.005 


0.37 ±0.03 


All 


-0.18 


0.03 


0.015 


0.000 


0.005 


0.34 ±0.03 


Si i 


-0.19 


0.02 


0.025 


0.015 


0.005 


0.33 ± 0.04 


Ki 


-0.32 


0.05 


0.045 


0.005 


0.005 


0.20 ±0.07 


Cal 


-0.41 


0.02 


0.035 


0.005 


0.005 


0.11 ± 0.04 


Sci 


-0.37 


0.06 


0.045 


0.000 


0.000 


0.15 ± 0.08 


Sen 


-0.29 


0.03 


0.005 


0.025 


0.015 


0.23 ± 0.04 


Til 


-0.25 


0.02 


0.050 


0.000 


0.000 


0.27 ±0.05 


Tin 


-0.31 


0.03 


0.005 


0.025 


0.010 


0.21 ± 0.04 


Vi 


-0.32 


0.02 


0.050 


0.005 


0.005 


0.20 ±0.05 


Cri 


-0.57 


0.03 


0.025 


0.000 


0.000 


-0.05 ± 0.04 


Mr I 


-0.73 


0.03 


0.020 


0.005 


0.005 


-0.21 ± 0.04 


Co I 


-0.43 


0.03 


0.000 


0.015 


0.010 


0.09 ± 0.04 


Nil 


-0.46 


0.02 


0.005 


0.015 


0.005 


0.06 ± 0.03 


Znl 


-0.30 


0.05 


0.030 


0.025 


0.005 


0.22 ± 0.06 



proximately true given the procedures used to constrain 
them. 

For Sc and Ti, spectral lines from two species, namely 
neutral and singly ionized, were available for our abun- 
dance analysis. In principle in both cases ionization bal- 
ance is satisfied within the lc uncertainties. The aver- 
age abundance for the two Ti species leads to [Ti/Fe] = 
0.23 ± 0.04. Nevertheless, it is clear that the line-to-line 
scatter for the Sci abundance is much larger than that 
for Sell. It is known that, for solar-type stars, Sci abun- 
dances are severe ly affected by non-LTE effects (e.g., 
iZhang et al.l2008T) whereas those obtained from Sc II lines 
are more reliable. The smaller error of the latter suggests 
that this could be the case also for giant stars. There- 
fore, for the Sc abundance we prefer to use that inferred 
from Sc II lines exclusively. 

Fig. [T^] shows our abundance results compared to 
Galactic chemical evolution trends based o n the analysis 
of large sam ples o f nearby dwarf s tars fr om lReddv et al.l 
(|2003T [2006h and IRamfrez et al.l (p007l) . These trends 
separate stars that are members of the so-call ed Galactic 
thin disk from those of t he thick disk (e.g . . iFuhrmannl 
[19981 iBensbv et al.1 12003 1 : IReddv et all 12001 l200l . The 
latter tend to have hotter kin ematics compared to the 
thin disk (e.g., iSoubiranl I1993D and based on their ob- 
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served properties (kinematics, chemical composition, age 
distribution, etc.), they are thought to be a different stel- 
lar population, separate from the thin disk, although 
its precise origin remains unknown. The (heliocentric) 
Galacti c space velocity comp onents of Arcturus, as de- 
rived by EamirezeEai] (JloTJ), are U = 25.2il.3kms" 1 , 
V = -119.0il.0kms- 1 , W = - 2.7 ±3.5 km s" 1 Usin g 
the membership formulation by iRamfrez et al.l ([2007) , 
we find that Arcturus has a probability of about 94 % of 
being a thick-disk member. 

As shown in Fig. [12l the chemical composition of Arc- 
turus is very typical of that of a thick-disk star in the 
solar neighborhood. The enhanced abundances we ob- 
tain for C, O, Mg, Al, Si, Sc, Ti, and V are also observed 
in ordinary thick-disk dwarf stars. The Ca abundance 
appears slightly low compared to the mean thick disk 
trend but it is still fully consistent with it considering the 
la star to star scatter. We note that our K abundance 
appears somewhat low compared to the m ean thin-disk 
trend (no K abundances were derived by iReddv et al.1 
120061 for thick-disk stars) , but we recall that the abun- 
dance of this element is based on the analysis of only one 
not so reliable feature. Relatively small abundance off- 
sets between Arcturus and the mean thick-disk trends in 
this comparison are not unreasonable given the impor- 
tant differences in stellar parameters between Arcturus 
and dwarf stars. 
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4. CONCLUSIONS 

Atmospheric parameters and elemental abundances of 
Arcturus, a primer for studies of red giant stars, are de- 
rived using high quality data and methods that mini- 
mize model uncertainties within the 1D-LTE approach. 
Limitations of these techniques are revealed in the incon- 
sistency between the iron abundance inferred from Fe I 
and Fe n lines and possibly also in the large error of Sc I 
abundances. These cases should be addressed using more 
sophisticated models, such as those including the impact 
of surface inhomogeneities and non-LTE. 

The elemental abundance pattern of Arcturus is typi- 
cal of that observed in ordinary nearby thick-disk stars 
and consistent with its sp ace velocity and relatively old 
age. iNavarro et al.1 (|2004l ) suggested an extragalactic ori- 
gin for the Arcturus group based on its angular momen- 
tum. Elemental abundance studies of stars in the Milky 
Way's satellite galaxies, however, show that extragalac- 
tic stars could be disentangled from the Galactic disk 
and halo stars using their detailed abundance patterns. 
For example, their [Mg/Fe] and [Ca/Fe] abundance ratios 
should be significantly lo wer than those observed in the 
solar neighborhood (e.g. JVenn et al.|[2004j : iTolstov et al.l 
2009). The elemental abundances we derive suggest that 
Arcturus was born in the Milky Way. 
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